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Abstract 

We consider an inextensible, semiflexible polymer or worm-like chain, with persistence length P 
and contour length L, fluctuating in a cylindrical channel of diameter D. In the regime D <C P -C L, 
corresponding to a long, tightly confined polymer, the average length of the channel (R\\) occupied 
by the polymer and the mean square deviation from the average vary as (R\\) = [l — a a (D / P) 2j/3 ^ L 
and {AR 2 } = (3 (D 2 / 'P)L, respectively, where a and /3 are dimensionless amplitudes. In earlier 
work we determined a and the analogous amplitude an for a channel with a rectangular cross 
section from simulations of very long chains. In this paper we estimate (3 and j3n from the 
simulations. The estimates are compared with exact analytical results for a semiflexible polymer 
confined in the transverse direction by a parabolic potential instead of a channel and with a recent 
experiment. For the parabolic confining potential we also obtain a simple analytic result for the 
distribution of Ri\ or radial distribution function, which is asymptotically exact for large L and has 
the skewed shape seen experimentally. 
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I. INTRODUCTION 



The statistical properties of biological polymers fluctuating in nano- or micro-channels 
have been studied in several recent experiments ■ For biological polymers the persistence 
lengths are typically tens of nanometers or even larger. When the channel diameter is 
smaller than the persistence length, the stiffness of the polymer plays an important role. 
The polymer is stretched out in the channel with little backfolding, and the length of the 
channel occupied by the polymer is only slightly shorter than its contour length. 

Measurements of the end-to-end distance of the polymer in a channel and its fluctuations 
provide information on the persistence and contour lengths of the polymer. This is of 
interest in studies of DNA fragments, for example, where sorting fragments of different 
ength is desired, or in determining the change in bending rigidity upon binding of proteins 

alsi- 

In this paper we consider the simplest model for a confined biopolymer - an inextensible, 
semiflexible filament or worm-like chain with persistence length P and contour length L in 
a cylindrical channel of diameter D. Here D is an effective diameter, equal to twice the 
maximum transverse displacement of the polymer from the symmetry axis of the channel. 



For this system the distribution of the end-to-end distance or radial distribution function 

H Q pjith, the change, replaced by a parabolic confining 



has been calculated theoretically 



potential, and studied with simulations [111 . Il2| 

We will mainly consider the regime D ^ P ^ L, corresponding to a long, tightly confined 
polymer. In this regime the length of the channel R» occupied by the polymer is essentially 
the same as the end-to-end distance. As discussed below, the distribution of R\\ is Gaussian 
and is completely determined by the mean value (R\\) and the mean square deviation (AR%). 
These two quantities have simple scaling properties, summarized in the next paragraph. Our 
goal has been to determine the dimensionless proportionality constants in the scaling forms 
with good precision, so that one has unambiguous predictions for the worm-like chain that 
can be compared with experimental data and used, for example, to determine the persistence 
length. 

In the regime D ^ P ^ L, the free energy per unit length of confinement Af, the 
average length of the channel occupied by the polymer, and the variance or mean-square 
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deviation from the average are given by 

Af = A 



(R\\) 
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as follows from scaling arguments of Odijk [3, arid a detailed microscopic analysis 



m 



16j. For a channel with a rectangular cross section with edges D x and D y , 
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Here A a , a , (3 , A n , and are dimensionless universal amplitudes, which do not 
depend on P, D, D x , and D y . 

The best estimates of the amplitudes in Eqs. (DQ), (j^J), (jlj), and (J5J) to date are 



A n = 1.1036 , 1.1038 ±0.0006 
a n = 0.09137 ±0.00007, 



A = 2.3565 ±0.0004. (7) 
a Q = 0.1701 ±0.0001 . (8) 

The first entry for A n in Eq. (JTj) was obtained by Burkhardt Ijj] , by solving an integral equa- 
tion numerically, which arises in an exact analytical approach. The other estimates are from 
our simulations 16j of very long polymers, with contour lengths up to L fa 1000 (2P 1 / 3 /} 2 / 3 ) 



where A ~ p 1 / 3 /) 2 / 3 is the characteristic deflection length introduced by Odijk 13]. Other 
estimates from simulations, compatible with these values but with larger error bars, are 
given in Refs. 
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Ref. 



2l| 



201 ] . and related results for a helical polymer in a cylindrical channel in 



The paper is organized as follows. In Section II the underlying theoretical framework is 
reviewed, and new estimates from simulations, 

(3 n = 0.0048 ± 0.0001 , O = 0.0075 ± 0.0002 , (9) 

for the amplitudes in Eqs. d3J) and Eqs. (jSJ), obtained with same method as in Ref. are 
presented. 



In Section III and the Appendix we consider the mathematically more tractable problem 
of a polymer tightly confined in the transverse direction by a parabolic potential instead of 
a channel with hard walls. Exact analytic expressions for each of the quantities Af, Rh, and 
(ARu) are derived. We find that (ARu) is overestimated by about 30 % if the potential 
parameters are chosen to reproduce L — (R\\ ) for a channel with hard walls. For the parabolic 
confining potential we also obtain a simple analytic result for the distribution of R\\ or radial 
distribution function, which is asymptotically exact for large L and for moderately large L 
has the skewed shape seen experimentally. 

In Section IV our predictions are compared with experimental results of Koster and 
Pfohl [4] for the radial distribution function of actin filaments in micro-channels. Section V 
contains closing remarks. 

II. THEORETICAL FRAMEWORK 

In the regime D ^ P ^ L, the line or filament by which we model the polymer is 
almost straight, without backfolding. Each such polymer configuration corresponds to a 
single valued function r(t), where (x,y,t) are Cartesian coordinates, and f— (x,y) specifies 
the transverse displacement of the polymer from the symmetry axis or t axis of the channel. 
Since the slope v = df/dt with respect to the t axis satisfies \v\ <C 1, the relation L = 
J fl " dt [1 + v(t) 2 ]^ 2 between the contour length L and the longitudinal length R\\ may be 



Here the two terms in square brackets are the bending energy per unit length and the 
confining potential per unit length, both divided by ksT. For a polymer in a channel with 
hard walls, V(r) takes the values and oo for r inside and outside the channel, respectively. 
According to Eq. ffTO]) . the average length of tube occupied by the polymer and its 



replaced by 



and the Hamiltonian 7-L of the worm-like chain 22| simplifies to 




(10) 




(11) 
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variance or mean square deviation are given by 

(R ll ) = L- 1 -£dt(v(t) 2 }, (12) 
(AJ2f ) = \ £ dt x £ dt 2 [{vih) 2 v(t 2 ) 2 ) - (v(h) 2 ) (v(t 2 ) 2 >] , (13) 

where AR\\ = R\\ - (R\\). 

For a tightly confined polymer in a channel with a rectangular cross section, the displace- 
ments of the polymer in the x and y directions are statistically independent. The partition 
function Z factors into a product of two partition functions Z x Z y , which only involve dis- 
placements in the x and y directions, respectively. This is a consequence of the additive 
property (d 2 f/dt 2 ) 2 = (d 2 x/dt 2 ) 2 + (d 2 y/dt 2 ) 2 in the Hamiltonian ffTTT) and the rectangular 
boundary, which does not break the statistical independence in the two transverse direc- 
tions. From this and from rescaling lengths according to x' = D~ l x, t' = (2P)~ 1 / 3 D X 2 ^ 3 1, 
it follows that the statistical averages on the right-hand sides of Eqs. ( fl2l) and (JTBl can all 
be determined from simulations of a long polymer with persistence length P' = \ confined 
to the two dimensional strip < x' < 1 in the (x',f) plane, as carried out in Ref. 16]. 

The statistical averages in Eqs. (|T2|) and (|T3|) can be expressed in terms of the variable 

C=^V<(0 2 , (14) 

where v' 2 = (dx' /dt') 2 . According to the scaling transformations in the preceding paragraph, 
/9P\ 2//3 1 r L /2P\ 2//3 1 f L 

<0= U) -lL 2 >=u) lL dt{v ' {ty2) - (15) 

As discussed in the final paragraph of the Appendix, the quantity £ defined in Eq. fTH|) 
is expected to follow a Gaussian distribution for sufficiently large L, with the mean value in 
Eq. f|T5l) and with variance w 2 given by 

ft' r t' 

t' 2 

The distribution determined from our simulations of polymers with values of t' up to 300, 
shown in Fig. 1, is indeed very nearly Gaussian, and the variance w 2 , as shown in Fig. 2, is 
in excellent agreement with the scaling behavior w 2 t' — > k for large t', where A; is a constant, 
expected [23J from Eq. (J6]). Substituting this relation in Eq. ( TT6|) and expressing the scaled 



w 2 = <(C " (C» 2 > = ^ / dt[ I dt' 2 [{v'M? <(t' 2 ) 2 ) - (v' x (t[) 2 } (v' x (t' 2 ) 2 }} • (16) 
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lengths in terms of the original variables gives 



2P 1 



Dl L 
2P 1 



L r L 



k = — - / dU dt 2 [{^(ti) 2 v x (t 2 ) 2 > - Mti) 2 )(v x (h) 2 )] 



JO 
L r L 



dt x / dt 2 [(^(ti)X(t 2 ) 2 }-(^(ti) 2 )K(*2) 2 }] • (17) 



£ 2 ^o .0 

According to our earlier paper [16], (C) = 0.2901 ± 0.0003, and from the data shown in 
Fig. 2 of this paper, we estimate k = 0.0382 ± 0.0010. Inserting these values in the relations 
a n = 2~ 5 / 3 (C) and /3 a = § k, which follow from Eqs. ©, ©, and (fl2>(fl?]) . we obtain the 
predictions for an and /?□ in Eqs. (jHJ) and (El). 

The entries for a and /3 Q in Eqs. (jHJ) and (j9j) follow, in a very similar way, from the result 
(Co) = 0.5400 ± 0.0004, where 

1 * 



obtained in Ref. 



Co = - / dt'v'{t')\ (18) 







16| from simulations of a polymer with longitudinal length t' and persis- 
tence length P' = | in a channel with a circular cross section of diameter D' — 1, and from 
the corresponding estimate fc = 0.0602 ± 0.0020, where = ((( — (Co)) 2 ) — > k Q t'~ l for 
large t'. 

III. POLYMER CONFINED BY PARABOLIC POTENTIAL 

Next we consider a polymer tightly confined in the transverse direction by a parabolic 
potential of the form 

V(f)= 1 -{b x x 2 + b y y 2 ) (19) 

instead of a channel with hard walls. The partition function Z(r,v;r\,,v* ;t) corresponding 
to the Hamiltonian (TTTTl with the parabolic potential energy ( TT9l) was evaluated for arbitrary 
values of the position and slope, (r,v) and (fo,Vo), at the polymer endpoints and arbitrary 



longitudinal length t in Ref. 



24J. 



The case b x — by of equ al potential parameters has been studied by Levi and Mecke 



and Thiiroff et al. llj, who calculated the distribution of R\\ or radial distribution 
function and compared their predictions with the experiments of Ref. j^, 4]. In this paper 
we consider distinct values of b x and b y , as is appropriate for rectangular channels with 
D x 7^ D y , and concentrate mainly on the large-L limit and on the prediction of the 6 
dimensionless amplitudes A , . . . , /3 a in Eqs. ([I])-©. 



Since the thermal averages in Eqs. (|T2|) and (|T3|) are integrated over the entire length 
of the polymer, the particular boundary conditions at the endpoints of the polymer are 
unimportant in the large- L limit. Straightforward calculations, given in the Appendix, lead 
to the results 

{k B TY l Af = 2- 1 ' 2 p- l ' i (b l J A + b l J A ) , (20) 
<J2||> =[i-2-^p-^ (6-V4 + b?*)]L, (21) 
(Ai2f ) = 2" 9 / 2 P" 5 / 4 (b^ 4 + b~^) L . (22) 

To obtain an approximate formula for the amplitude fin for a channels with hard walls and 
a rectangular cross section, defined in Eq. (jSJ), we choose the parabolic potential parameters 
b x and b y in Eq. (j2ip so that the average longitudinal length in the channel, given by 
Eq. (jSJ), is reproduced, term by term. Substituting these potential parameters in Eq. ( 122|) 
and comparing with Eq. ([6]) leads to a formula for fi n in terms of This calculation and 
a similar one for the channel with a circular cross section lead to the relations 

fin = 8a D , fi = 2al . (23) 

We note that Eq. ( 12 3 p also follows from choosing the parabolic potential parameters in Eq. 
( 122|) to reproduce {AR^} in Eqs. (E} or (J6]), substituting these potential parameters in Eq. 
f l2T|) . and comparing the result with Eqs. (J2]) or (jSJ). 

Substituting the values of a a and a in Eq. (JSJ) into Eq. fT23|) . we obtain the predictions 
/?□ = 0.00610 ± 0.00002 and fi = 0.00984 ± 0.00002, which are 27 % and 31 % larger, 
respectively, than our estimates from simulations. Thus, we see that calculations in 
which the hard wall potential of is replaced by a softer, parabolic confining potential tend 
to overestimate the endpoint fluctuations (AR 2 ) if the potential parameters are chosen to 
reproduce L — (R») for a channel with hard walls. Similarly, if the potential parameters 
are chosen to reproduce (AR 2 ) for a channel with hard walls, the quantity L — (R\\) is 
underestimated. 

The asymptotic forms of both L — (R\\) and (AR„) for a polymer in a channel with hard 
walls are correctly reproduced if not only the potential parameters, but also the persistence 
length P of the equivalent parabolically confined polymer is properly chosen. Setting Eqs. 
f l2T]) and fT22|) . with P in place of P, equal to the corresponding expressions ©, ©, (JSJ), 
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and ©, and solving for P, we obtain 

p = I^P^ p=Ap (24) 

for the rectangular and circular channel cross sections, respectively, where the same combi- 
nations of exponents occur as in Eq. fl23|) . Substituting the values of a a , a , 0n, and (3 Q 
from Eqs. (jHJ) and (EJ) in Eq. (1241) . we find that the persistence length P of the equivalent 
parabolically confined polymer is 21 % and 24 % smaller than the persistence length P of 
the polymer in the rectangular and circular channel, respectively. 

Finally, in the Appendix we derive simple analytic results, in terms of "inverse Gaussian" 
functions, for the radial distribution function of a polymer confined by a parabolic potential 
in the moderate to large L regime. The predictions, given in Eqs. (lA14j) - (1A16[) . OA20p . and 
(lA2ip . with £ — L — R\\, are compared with experimental results for polymers in channels 
in the next section. 



IV. COMPARISON WITH EXPERIMENT 

Experiments on unconfined filaments of the biopolymer actin (see Ref. and references 
therein) have yielded estimates of 8 to 25 /im for the persistence length. With fluorescence 
microscopy Koster, Pfohl, and coworkers have measured the radial distribution of 

actin filaments with contour length L = 21 fim in channels with rectangular cross sections 
with depth D x = 1.4 /zm and widths D y = 1.5, 4.0, 5.8, and 9.8 /im. Comparing their 
experimental results for the radial distribution function, shown below in Figs. 3 and 4, with 

n 

the theoretical prediction of Levi and Mecke [10[ for a parabolic confining potential, Koster, 
and Pfohl [4] find good agreement, for all four channels, with the value P = 13 /im. 

Since L is only moderately larger than P, the above experimental parameters do not 
clearly satisfy D x ,D y <C P <C L, the condition under which our predictions for (R\\) and 
(Ai?i?) apply. Nevertheless it is interesting to compare the experiments with our predictions 
for the scaling regime. 

As discussed above and in the last paragraph of the Appendix, the distribution of Ru is 
expected to be Gaussian in the scaling regime, with mean value and variance given by Eqs. 
©, (J6]), (|8]), and Qj. Using these relations and the above experimental values of L, D x , and 
D y to determine the mean and variance as a function of P, we have carried out least square 
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fits of the experimental results to Gaussian distributions for all four channels, varying P to 
optimize the fits. This leads to the results shown in Fig. 3, and the estimates P = 7.61, 
11.1, 14.1, and 10.1 //m for the the channels with widths D y = 1.5, 4.0, 5.8 and 9.8 /xm. 
The first two of these estimates are expected to be the most reliable, since the condition 
D x , D y P <C L is more nearly satisfied. 

We have also carried out fits of the experimental results in which both (R\\) and (Ai2|) 
are treated as fit parameters. In the large- L limit these quantities yield two independent 
predictions, 

/ ^2/3 2/3\ 3 / 2 

D 2 + D 2 

p=k ji^iL' (26) 

for the persistence length, which follow from solving Eqs. (JSJ) and fl6]) for P. 

Least square fits of the same experimental data to the inverse Gaussian distribution, 
given by Eqs. (IA14j) and (IA18I) . with both the mean (R\\) and variance (ARn) adjusted 



to optimize the fit, are shown in Fig. 4. Of course, the two-parameter fit reproduces the 
experimental distribution more closely than the one-parameter fit in Fig. 3. Both the inverse 
Gaussian distribution and a convolution of inverse Gaussian functions, as described in the 
Appendix, have the skewed form seen in the experimental data and lead to nearly the same 
results. 

The fits shown in Fig. 4 lead to the estimates P = (7.0, 2.8), (9.5, 3.6), (10.8, 4.1), and 
(7.2, 3.2) in /xm for the channel widths D y = 1.5, 4.0, 5.8, and 9.8 jum, where the first and 
second numbers in parenthesis follow from substituting the mean and variance from the best 
fit in Eqs. ( )25l) and ( |26|) . respectively, with an and /?□ given by Eqs. (IE]) and (J9j). All of 
these estimates are smaller than the values P = 13 fim and P = 15 ± 3 /mi proposed in 



Refs. {4, [lOj, respectively, and for each channel the estimate based on Eq. ( 126]) is only 3 
or 4 /im, less than half of the corresponding estimates based on Eq. fl25|) . Determining the 
mean and variance by fitting the experimental data to an ordinary Gaussian distribution 
instead of an inverse Gaussian distribution or by evaluating the mean and variance directly 
from the experimental histograms without assuming a particular distribution leads to quite 
similar estimates. 

Finite-size corrections probably account, at least in part, for the discrepancy in the esti- 



mates of P based on Eqs. (125]) and (|26|) . with the smaller estimate coming from Eq. (126]) . 
As the contour length L increases and the polymer is tightly confined over a greater fraction 
of its length, (Ai?|)/L approaches its limiting value from above, so that P, as given by Eq. 
(T2"6"]) . approaches its limiting value from below. In Fig. 2 the lower and upper curves in 
the inset show the finite size corrections for polymers with one free end and two free ends, 
respectively, with the latter case corresponding to the experiment. For P = 10 /im, D x = 4 
/mi, L = 21 /im, the rescaled length t' = (2P)- 1 ^Dl /3 L is about 3.1, and for this value 
of t', (ARn)/L oc w 2 t' is seen to be about 50 % larger than its large t' limit. The actual 
finite-size corrections are expected to be even larger than this, since Fig. 2 is based on the 
Hamiltonian f llip . which is equivalent to the worm-like chain for small slopes \dr/dt\ <C 1 , 



22]. 



but for larger slopes overestimates the bending energy 

In comparing the estimates of P from Eqs. (125]) and ( 126]) . one should keep in mind that 
the prediction of Eq. ( 125]) is extremely sensitive to the experimental uncertainty in the 
normalized mean (AR\\)/L, since this quantity is close to unity for a long tightly-confined 
polymer, so that the denominator in Eq. (125]) nearly vanishes. For example, increasing 
(AR\\)/L from the value 0.93 by 3 % more than doubles the estimate of P. In view of this, 
the disagreement of the numerical estimates based on Eqs. (125]) and f]26]) mentioned a few 
paragraphs above is not so surprising. One advantage of Eq. ( 126]) over Eq. ( 125]) is that the 
relative uncertainties in P and (ARu)/L are the same. 



V. CONCLUDING REMARKS 



In Ref. [16|] and this paper we have determined the universal amplitudes a , (3 , and 
/?□ in the scaling forms ([2]), ([3]), ([5]), and ([6]) for the worm-like chain in cylindrical channels 
with good precision from simulations. We hope the results will be useful in analyzing ex- 
periments. Combining measurements of (R\\) and (AR„) and our predictions, one obtains 
two independent predictions for the persistence length P, which can be checked for consis- 
tency. We recall that (ARn) may be determined by measuring the isothermal extension of 



251 ] as well as by direct observation of 



a polymer in a channel placed under a weak tension 
the endpoint fluctuations. 

We have also derived exact analytic results for a polymer confined by a parabolic potential 
rather than a hard wall and shown that (AR2) is overestimated by about 30 % if the potential 



10 



parameters are chosen to reproduce L — (R\\) for a channel with hard walls. 

Finally, we have compared our predictions for the scaling regime with the experimental 
data of Ref. 4| for the radial distribution function. The comparison points to a per- 
sistence length smaller than the values 13 fim and 15 ± 3 fim reported in Refs. |4| and 
[lo| . respectively, but the experimental parameters are at the edge or outside the scaling 
regime, and significant corrections to scaling are expected. For a more conclusive comparison 
with our results, we would welcome experiments that probe deeper into the scaling regime 

D x , Dy^P^L. 



Appendix A: Calculational details for parabolic confining potential 



For the Hamiltonian (TITl) with the potential energy ffl9l) . the polymer partition function 



for a polymer in the three dimensional space (x, y, t) factors in the form 



Z ^ (clt,, br,', o 1( , b qi 



zf'(a x ,b x )Z 



(2), 



a y , by) 



(Al) 



Here 



r(2), 



a, b) 



.Do; exp 



dt 



P 



d 2 x 
~d& 



+ a(^ + bx 2 



(A2) 



is the partition function of a worm-like chain in two spatial dimensions (x, t), with a parabolic 
confining potential. 

In Eqs. flAlj) and flA2j) . auxiliary fields a x and a y have been introduced for conveniently 
generating correlations of dt v\ and dt Vy by differentiation. The auxiliary fields have 
a physical interpretation related to tension. If one end of the polymer is fixed and the other 
end is free to move but subject to a force or tension r applied in the longitudinal direction, 
the corresponding potential energy —t{R\\ — L) ~ | dt v(t) 2 , where we have used Eq. 
(flOl) . is included in the Hamiltonian and contributes to the Boltzmann factor. Comparing 
with the partition functions in Eqs. ( lAlj) and ( 1A2[) . we see that a x = a y = r/kBT. 



For calculating "bulk" properties of long polymers that are independent of the detailed 
boundary conditions at the ends, the periodic boundary condition r(t) = r(t+L) is especially 
convenient. With the substitution x(t) = L~ x l 2 ^2 q x q e iqt , Eq. ( 1A2I) takes the form 



'(2), 



a, b) 



Dx Y\ exp 

9 



-- (Pq A + aq 2 + b) x q 



X- 



(A3) 
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The subtracted free energy Af^(a,b), defined by 

Af^(a,b) 



L' 1 In 



Zf{a : b)/Zf (0,0) 



(A4) 



may be evaluated by standard Gaussian integration techniques 26j and is given by 



Af( 2 \a,b) 



[°° dq Pq* + aq 2 + b = /2 /4 /4 / 



1/2 



(A5) 



Jo 2vr Pg 4 V 2V&P, 

The right-most expression in Eq. ( 1A5I) also follows readily from the path-integral ap- 
proach of Ref. J24I], according to which the partition function of the polymer with fixed 
endpoints and endslopes has the expansion 



Z i - 2) (x,v : x ,v \L) = ^27p v (x,v)ip v (x , -v )e~ 



E v L 



(A6) 



analogous to a quantum mechanical propagator. The eigenvalues and eigenfunctions are 
solutions of the L-independent Fokker-Planck equation 

d 



'^~2P~J^ + \ bx2 + \ av2 )M X > V ) = Ei/j(x,v) . 



(A7) 



The dominant contribution for large L in Eq. (IA6j) comes from the ground state, which has 
eigenfunction ip (x,v) and eigenvalue E , where E = (k B T)~ 1 Af( 2 \a,b), as follows from 
Eqs. ( ]A4j) and fl A6[) . According to Ref. 24| . ^0(2^) has the Gaussian form ip (x,v) = 
Aexp (—Bx 2 — Cxv — Dv 2 ). Requiring that this expression satisfy Eq. ( 1A7I) determines Eq 
and the constants B, C, D, and Eq. ( 1A7I) . yielding 



= Aexp [-(6P) 1/2 P a; 2 + (oP) 1/2 xi; - PP v 2 ] 



(A8) 



with Po given by the right-most expression in Eq. ( 1A5I) . 

Setting a=0 in Eq. (1A5j) and including the contributions from displacements in both the 
x and directions into account, we obtain the free energy per unit length of confinement in 
Eq. ([20]), which is consistent with Eq. (16) of Ref. [24]. 

From Eqs. (1A"2]) . (TAI) . and (fA5l) . 



9 A/( 2 )(a,6) 1 



-V / dt (v(t) 2 } = 2 - 5 / 2 r 1 / 4 p- 3 / 4 f 1 + 



-1/2 



(A9) 



9a fc B T 2LJ Q V 2V6P, 

To calculate the average longitudinal extension (R\\), we set a=0 in Eq. (1A9[) . substitute the 
result in Eq. ffl2|) . and include the contributions from transverse displacements in both the 



12 



x and y directions. This yields the expression for the average longitudinal extension given 
in Eq. (j2T|). 

Similarly, from Eqs. (TM]), O, and (TA5|l . 



d_Y AfW(a,b) _ 1 
<9a 



4L 



2 -9/2 6 -3/4p-5/4 1 + 





-3/2 



(A10) 



To obtain (Ai?,? ) , we set a=0 in Eq. (lAlOj) . substitute the result in Eq. ( TT3|) . and include 
the contributions from tranverse displacements in both the x and directions. This yields 
the expression for the average longitudinal extension given in Eq. ( 1221) . 
It is straightforward to derive the complete distribution function 

-L 



v(ty 



(All) 



from which the above moments follow. Its Laplace transform is given by 

'■L 



P^ (s; a, b) = J™ di e -*P(£; a, b) = /exp f-~ ^ 



(A12) 



where the average is to be carried out with the same Boltzmann weight as in Eq. (1A2I) . 
Thus, 



P (2 \s-a,b) 



Z?\a + s,b) 



Z?\a,b) 



exp 



—L 



fW(a + s,b)-f 2 \a,b) 



k B T 



(A13) 



where we have made use of the definition flA4p . Substituting Eq. fl A5j) in Eq. (IA13P and 
equating the inverse Laplace transform, we find that P« K lM) is given by the "inverse 
Gaussian" or Wald distribution I2j 



Pi, 



,(0 



1 



<0\ 

/ 2^(AO U/ 



3/2 



exp 



(0 (£ - 

£ 2<AO J 



< £ < oo 



(A14) 



where 



(0a, b = 2-^b'^P-^ (l + 

a 



-1/2 



<Af) 0)6 = 2-^b^P-^ (l + 



-3/2 



L 



(A15) 
(A16) 



are the mean and variance of the distribution, respectively, consistent with Eqs. (1A9|) and 
f lATOj) . Note that inverse Gaussian distribution vanishes as £ approaches zero, as expected 
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from Eq. (lAllj) . reflecting the fact that the end-to-end distance R» of the polymer cannot 
exceed the contour length. 

Since the mean and variance in Eqs. flA150 and flA160 are both proportional to L, the 
inverse Gaussian distribution (IA14D reduces to the ordinary Gaussian form 



PgaUSs(0 ~ v^(AO 6XP 



2(Ae) 



oo < £ < oo , (A17) 



in the large-L limit. 

The above results for a polymer in a two dimensional space (x, t) are easily generalized 
to three spatial dimensions. In Eq. flAHI) . the quantity v 2 is replaced by v^+v 2 , so that 

£, = L — R\\ , (A18) 

in agreement with Eq. (jTUJ), and Eq. flA13[) is replaced by 

P^(s-a b-a b ) = + g; hx) Z ^ ){flv + hy) (A19) 

Z K L J (a x ,b x ) Zl J (a y ,b y ) 

Accordingly, the inverse Laplace transform is given by the convolution 

^ 3) (0 = fd?p( 2 \Z-Z>;a x ,b x )pW(Z';a y ,b y ) , (A20) 
Jo 

where each of the factors P^ in the integrand has the inverse Gaussian form (1A14I) . with 
mean and variance defined by Eqs. (1A15j) and (1A16I) . 

In the case of cylindrically symmetric potential parameters 
appropriate for a channel with a circular or square cross section, the convolution in Eq. 
( 1A20I) can be evaluated (or circumvented). The corresponding distribution also has the 
inverse Gaussian form 

^ (3) (0 = iW»(0> (A0 = 2(Ae) a ,b, (A£ 2 > = 2(Ae) a ,b , 0<e<oo, (A21) 

in terms of the distribution (1A14I) and the mean and variance defined in Eqs. ( 1A15I) and 
(IA16|) . 



In the large-L limit, in which p( 2 )(£;a, b) becomes Gaussian, the distribution functions 
Po\Cj and Pj 3 ^(£) both take the Gaussian form (1A17|) . with mean (£) = {£,}a x ,b x + {Oa y ,b y 
and variance (A£ 2 ) = (A^ 2 ) axt b x + (AC i 2 ) ayt b y defined in Eqs. ( 1A15I) and ( 1A16I) . as is consistent 
with Eqs. (ED) and (|22l). 
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Like Eqs. (Ell) and flAH]), our predictions flA20|) and flMTj) for the distributions P n (f) 
and P (£) in terms of inverse Gaussian functions are really only exact in the large- L limit, 
in which the ground-state contribution to the sum in Eq. (IA6[) dominates. However, for 
moderately large L the distributions also work quite well, reproducing the skewed form of 
the radial distribution observed experimentally and calculated theoretically in Refs. lOL 111]. 
This is shown in Section IV, where our results are compared with recent experimental data 
of Koster and Pfohl [4] for the radial distribution function. 

Finally we argue that the distribution of Ru becomes Gaussian in the large-L limit not 
just for the parabolic potential, but for general confining potentials, including the hard-wall 
potential. To see this, note that for a general confining potential, the Laplace transform of 
the distribution function, defined as in Eqs. flAlip and (IA12[) is related to the free energy 
per unit length f{a) by 

fW(a + s)-fW(a) 



p(2) ( 



s; a) 



exp 
exp 



-L 



k B T 



1 



-(i)s + -(t\e)s 2 + LO{s*) 



(A22) 



analogous to Eq. (IA13I) . Here we have expanded f{s + a) to second order in a, relating the 
expansion coefficients to moments of £, as above. With the substitution s = iy the inverse 
Laplace transform of Eq. flA13|) takes the form 



1 f 00 

P {2 \^a) = — J dyexp 



z(£-(£>h--<A£ 2 )y 2 + LO((^) 3 ) 



(A23) 



Treating the O {{iy) 3 ) term in square brackets perturbatively, one finds a negligible contri- 
bution, for large L, to the Gaussian distribution ( 1A17I) implied by the first two terms. 
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FIG. 1: Distribution of the quantity £ = t' 1 J * dt' for a rescaled polymer with persistence 
length P' = \ and longitudinal length t' on a two-dimensional strip of width 1. The curves 
correspond, from bottom to top, to t' = 100, 225, 400, 625, and 900. 
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FIG. 2: Dependence of w 2 t' on t'~ , where w 2 is the variance or mean square deviation of the 
distribution in Fig. 1. The straight line shows the best two-parameter fit of the data from t' = 11 
to 300 to the functional form w 2 t' = k + it'" 1 , for which k = 0.0382 and I = 0.0129. The round 
points are simulation results for a polymer with one end fixed in the middle of the strip with slope 
vo = and with the other end free to fluctuate. The square points are results for a polymer with 
both ends free to fluctuate, as in the experiments. As seen in the inset, the finite size corrections 
to the limiting value for large t' are greater in the case of two free ends. 
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? IG. 3: One-parameter least squares fit of the experimental results (histogram) of Koster and Pfohl 
J] for the radial distribution function of actin filaments in channels with cross section D x x D y 
with the expected Gaussian distribution (dashed curves) for D x ,D y CPcL, The histograms 
and dotted curves are normalized to unit area. The mean and variance of the Gaussian curves were 
determined from Eqs. © and ©, using the estimates of an and /?□ in Eqs. flSJ) and ©, for the 
experimental parameters L = 21 fj,m, D x = 1.4 /im, and -D^ =1.5, 4.0, 5.8, and 9.8 /im. Choosing 
the persistence length P to optimize the fit, yields the estimates P = 7.61, 11.1, 14.1, and 10.1 /im 
for the channels of width 1.5, 4.0, 5.8, and 9.8 /im. 



20 





40- 



_Q 



10- 



1.4x5.8//m 2 









0.6 



0.7 



0.8 
fl||/L 



0.9 



1.0 




FIG. 4: Two parameter fit of the same experimental data for the radial distribution function as in 
Fig. 3 to the inverse Gaussian distribution, given by Eqs. (lATIj) and (jXlSj) . with both the mean 
(R\\) and variance (Afly ) chosen to optimize the fit. This leads to the estimates P = (7.0, 2.8), 
(9.5, 3.6), (10.8, 4.1), and (7.2, 3.2) in /iin for the channel widths D y = 1.5, 4.0, 5.8, and 9.8 /xm, 
where the first and second numbers in parenthesis follow from Eqs. (|25p and ([26]), respectively, 
with an and f3n given by Eqs. ([S]) and Q. 
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